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Pattern recognition of MorayeS HEREC fluctuations, representing 
plumes, occurring in data recorded in the atmospheric boundary layer 
was performed utilizing predefined characteristics of these phenomena. 
An averaged, normalized plume was computed based on the character- 
istics of ten previously identified sihneae. This standard plume con- 
sisted of two components or curves and a component-by-component 
correlation was made of this standard plume with various data sets, 
each set containing eaeeeaiwe information recorded at two different 
heights. A product curve was then formed from the two correlation 
curves and used to identify plumes automatically. A plume was said 
to have been detected if the product curve exceeded a previously 
selected threshold. The accuracy of this detection method was de- 
termined by visual inspection of the same records by a meteorologist. 
Approximately 73 per cent of all plume occurrences in the data sets 
were detected which indicates the appropriateness of this method of 


pattern recognition for plume detection. 
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I INTRODUCTION 


Patterns may be broadly seemed into two major groups; man- 
made and natural. Man-made patterns, suchas the characters inan 
alphabet, consist of symbols which are designed and produced by 
humans. Natural patterns, suchas cloud formations, are those found 
in nature to which man assigns ae BASE 

Pattern recognition 1s an emerging field of study with many 
applications to both man-made and natural patterns. It is nota co- 
hesive discipline in its om right, however, buta vast collection of 
highly varied problems. Any technique which contributes to the 
solution of any of these problems can be considered an integral part 
of pattern recognition. These techniques may be statistical, adaptive, 
or heuristic depending upon the special characteristics of the pattern 
type under consideration. 

A natural pattern of interest occurs within traces of temperature 
data recorded inthe atmospheric boundary layer. This pattern ap- 
pears as an intermittent asymmetrical saw-tooth "pulse" of large 
amplitude and always occurs as a positive fluctuation in the ambient 
temperature. The Sore nee of these pulses in data of this kind has 
been reported by Taylor (1958), Priestly (1959) and Webb (1965). 
Taylor suggested that the pulses are due to organized thermal struc- 
tures, or "plumes", contained within the complex turbulent flow. A 
schematic of such a plume appears in Figure l. 

Temperature inhomogeneities in the non-ionized atmosphere, such 


as those associated with plumes, are among several mechani: ms which 
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distort the propagation of an electromagnetic wave over a line-of-sight 
path. Strohbehn (1968) placed these distortions into three groups: 

(1) absorption; (2) refraction; and (3) random scattering due to random 
fluctuations in the dielectric constant which causes variations in the 
amplitude, phase, angle-of-arrival, and polarization of the wave. 

The last of these distortions is of interest because the final 
accuracy of radar systems, communication networks, and other 
systems may be limited by propagation characteristics. Furthermore, 
propagation effects that are negligible at lower frequencies may be 
dominant in higher frequency systems such as in the millimeter and 
optical bands. 

Plumes are also of interest for their role in the exchange of heat, 
moisture, and momentum between the earth's surface and the atmo- 
Spnere. Pricstly (1959} contends that knowledge of the role of plumes 
in transferring heat and water vapor to the atmosphere is essential 
to properly define boundary fluxes of these quantities. 

Many scientific disciplines require statistics describing plume 
phenomena. The versatility and usefulness of pattern recognition 
techniques and results underscore its value in this data-gathering 


task. 








I epEsCRIPTION OF pose 


All measurements of temperature fluctuations used in this study 
were originally recorded by Kenneth L. Davidson in 1968 under a grant 
administered by the Office of Naval Research. Davidson (1970) describes 
in detail the acquisition and processing of these data. The techniques 
and procedures used by him in this data-gathering task were studied to 
insure that the appearance of plumes represented actual geo-physical 
phenomena and were not artifically introduced via equipment irregular- 


ities. An abridged version of his description is presented in this chapter. 


Pa Ele LD SITE 

All data were obtained on Lake Michigan from a fixed U. S. Lake 
Survey Research Tower located in 15 meters of water 1.6 km from 
shore near Muskegon, Michigan (Figure 2). The uncomplicated bottom 
and shoreline features present in this area were reasons for selecting 
this site. The lake is surrounded by first order weather stations which 
permitted the acquisition of complementary meteorological information 


during measurement periods. 


fe DESCRIPTION OF TOWER 

Figure 3 shows the general construction of the research tower and 
its instrumentation during acquisition of the data used in this study. 
The structure was designed to achieve stability and to minimize dis- 
turbance to the airflow and water motions. The various sensor mount- 
ings on the tower enabled collection of turbulent velocity, wind, and 


wave measurements as well as temperature data. Temperatic 2 


measurements in the air were made simultaneously at two heights. 
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G. EQUIPMENT SENSORS FOR TEMPERATURE MEASUREMENT 
A resistance thermometer was used to measure temperature 
fluctuations. The sensor, a 30 ohm, .00015 inch tungsten wire, .16 
inches long, is part of a constant-current circuit and with a small 
filament current (approximately 2 ma) serves as a resistance thermo- 
meter. The bridge unbalance caused by a temperature fluctuation, and 
hence a sensor resistance change, was amplified 2500 times for the 
output signal. A block diagram of the bridge circuit is illustrated in 


Figure 4. 


>. RECORDING SYSTEM 

All data we re recorded on magnetic tape in analog fees using 
frequency modulation. A seven channel recorder was used in this 
process, with all recordings made at a tape speed cf 3 3/4 inches per 
second. The corresponding recorded frequency was D.C. to 625 Hz. 
Two of the seven channels available were allocated for the two tem- 
perature sensors. A schematic of the complete recording system 
(recorders and preamplifiers) appears in Figure 5. 

Calibration signals were recorded on all channels at the beginning 
and end of each tape and at i point where recording was interrupted. 
These signals served as indicators-‘of recorder Ae as well as noise 
level and were used during digitizing procedures. Constant monitoring 
of amplified and pre-amplified signals was performed witha dual 
channel oscilloscope in order to evaluate the electronic outputs and to 
insure that the input signals would not exceed amplifier or recorder 


input limits. 
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currrent); Flow Corporation System, Series 900. 
(after Davidson, 1970) 
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Figure 5, Schematic of recording system components. 
(after Davidson, 1970) 
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He DATA REDUCTiGi 

Steps performed to digitize the data are summarized in Figure 6. 
Because the recorded signals could range only from -2 to +2 volts and 
the analog-to-digital converter hada range from -100 to +100 volts, a 
gain of 50 was applied to all signals. Calibration signals recorded 
before, during and after measurement periods were used to adjust for 
Bain and drift errors in the recorder electronics. 

Data from all channels were lowpass filtered through identical 
second-order Butterworth analog filters. Figure 7 shows the analog 
circuit and the energy response function. A sampling rate of 300 
points per second was used in digitizing. 

The temperature fluctuations of interest were those with frequencies 
below 10 Hz, so additional numerical filtering was performed. A feed- 
forward fiiter was ernpioyed and Figure 8 shows its magnitude - squared 
frequency response characteristics. Also included are the tap weights 
and transfer function, H(w), for the filter. Since the digital filter output 
was only computed for every tenth point in the original data, the data out 
of this numerical filter consisted of 30 samples per second ofa signal 
with frequency components ranging from D.C, to about 10 Hz. Since 
plumes appear to have no significant frequency components above 2.5 
Hz, the data were filtered again (numerically) to eee every sixth 
point resulting ina filter output of 5 points per second. This was the 
data rate used in this study. The response curve for this filter appears 


merigure 9. 
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Figure 6. Flow chart of procedures in analog to digital 
processing. (after Davidson, 1970) 
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employed in analog low-pass filtering. 
(after Davidson, 1970) 
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4) H(w) = h. + 22 h_cos(ndtw) 
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where ® = radian frequency; 


At = time interval betweeen data points 
41) h(n t) = pee ae sin(wendt) 


Ee 
T= (w.-0,)© (ndt)® 


1 to insure no change in mean 


where: an = 
= number of weights (21) 
w, = cut-off frequency (w fem = 10 Hz) 
w, = terminal frequency (w/e = 20 Hz) 
Figure 8. . Frequency response, (A/A vee curve and weight 


functions employed in numerical low-pass filter. 
(after Davidson, 1970) 
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iil. THEORETICAL CONSIDERASHIONS 


Visual examination of temperature traces containing plumes indicates 
a striking resemblance to displays of communication signals entrained 
in noise. Since the optimum detector of such signals in white Gaussian 
noise is a matched filter, attention was directed toward this technique 
as a detection method for plumes. Accordingly, the data used in this 
study were analyzed to determine their distributions. The results of 


this analysis are summarized in Table I. 





Data Height Standard ,Skewness: Kurtosis: 
Period Level (m) Devaation™ -3,,-2%82 a4 ,,.cM 
en em) e /(€°) 
j ra, Os Oe su IAM Oe | 
10, 0.16 0.38 Coe 
9 10 Ora 0436 2. 66 
rae ® O-arS Ont02 fers) 
3 18, 0.18 O22ie (ay Tire: 
fare, 0.16 0.49 cago 
‘i ares, Oa Ore 2 Ze ow 
4.-0 G25 -0.09 (aie, 





Nr 


*e is the average value of the fluctuations 


Table I. Statistics of Temperature Data 


The tabulated values of skewness in Table I indicate that the data 
do not possess a Gaussian distribution. The noise itself may be 


normally distributed, however, and the combination of plume plus 


fags 








noise may skew the distribution. Therefore, the use of the matched 
filter approach to plume detection was considered. 

ihe matched filtey coe applied to pattern recognition can be 
described as a form of template matching. The categorization process 
is then equivalent to a component-by-component correlation of an un- 
known sample with a known sample or pattern. Although correlation 
methods can be shown to be optimal only under certain very restrictive 
conditions on the distributions in question, these methods were con- 
sidered the logical point at which to initiate attempts at solution of the 
peeerae recognition problem. Furthermore, since the plume detection 
problem is a binary one, either the plume exists in the random turbu- 
lent flow or it does not, the correlation method is particularly ap- 
mprepriate. If the measurement space contained many distinct patterns 
to be categorized, then different approaches for delineating hyper- 
planes and establishing clustering criteria should be considered. 
Correlation techniques, although Snes nel simple, occupy a 


valid position in the framework of pattern recognition. 
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IV. EXPERIMENTAL PROCEDURE 


The correlation method of pattern recognition requires a known 
template sample to compare against unknown samples. Applied to 
plume detection, the template becomes an averaged, ''standard" plume. 
Salient plume characteristics considered in the computation of this 
standard plume are listed below: 

(1) a gradual rise of the temperature data followed by a 
sharp drop, producing an asymmetrical saw-tooth or 
triangular trace, 

(2) simultaneous occurrence of the pattern in temperature 
traces from two measurement levels arranged ina 
plumb line, 

(3) four to six seconds duration of each pattern, 

(4) separation of patterns by quiescent periods, 

(5) fluctuations contributing to the pattern are always 
positive with respect to the recorded trace's 


preterred mean value. 


CAL-COMP plots were obtained for one 18 minute data set and, 
using the above Sheer ristics, a total of ten plumes were identified. 
The exact point where the sharp drop occurred in each plume was 
determined and a computer program was developed to compute an 
averaged, normalized, standard plume. The flow chart of this pro- 
gram appears as Figure 10 and the standard, normalized plume 


defined by this procedure is shown in Bicwre oll ~ 
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PEUViTE 
of 40 plumes to mean QO, 
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Figure 10. Flow chart of program to compute standard plume and 
correlate it with data sets. 
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This standard plume, consisting of two traces to correspond with 
the two heights at which the temperature fluctuations were measured, 
was correlated with the data in four ways: 

(1) correlation of the first component of the standard 
plume with the top trace of the temperature data, 

(2) correlation of the second component of the standard 
plume with the bottom trace of the temperature data, 

(3) addition of both correlation traces, 


(4) multiplication of both correlation traces. 


The correlation traces obtained were positioned next to their respective 
temperature traces and the varying degrees of correlation magnitude 
were noted in relation to plume occurrence. 

Correlation plots were initially obtained on the same data set 
from which the standard plume was computed. ‘This data set can be 
considered as ''training set'’ data. The standard plume was then cor- 
related with various ''test set'' data. To establish confidence in the 
validity of the standard plume, two different test sets were used. The 
first consisted of data acquired the same day and at the same height 
levels as the training set. Meteorological conditions were also uni- 
form during both periods. The second test set data was obtained ona 
different day, at different heights, and with different meteorological 
conditions prevailing. | 

Kenneth L. Davidson, Assistant Professor of Meteorology at the 
Naval Postgraduate School, subjectively analyzed the traces obtained 
and identified plume occurrences. Thresholds for low, medium and 
high probability of detection were then determined from the correlation 


plots. 
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VV. PRESENTARION OF DATS 


A total of four measurement periods, each eighteen minutes in 
length, were considered in this study. The heights of the two tem- 
perature sensors above the mean water level were the same for 
periods 1, 2, and 3 but were different for period 4. Data from period 
] were used as the training set, data from periods 2 and 3 were com- 
bined to form the first test set and data from period 4 were used as 
the second test set. Pertinent information concerning these periods 


is summarized in Table II? 


Significant 
Time Measurement maxe Height 
Period Date (bes i) J.evel (cm) 
1/2 
(meters ) 
September 27 
: (1968) — 1120-1138 1.0 and 2.0 42 
September 27 
2 (1968) 1140-1158 1.0and2.0 46 
3 PORE” 2) loa I. © andl 2, C 44 
(1968) 
September 26 
4 (1968) Veo overage clectae 38 
oe 
Hy /3 = four times the standard deviation of the water 


surface. The value of Hj/;3,is a good estimate of 
the average height (crest 143 trough) of the highest 
one-third waves inthe record. 


Table IJ. Description of Sampling Intervals and Physical 
Conditions During Data Acquisition. 
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A. WEATHER CONDITIONS DURING SAMPLING INTERVALS 

Light to moderate wind speeds prevailed during the first three 
periods. The mean coe velocity at the 4 meter level at the onset of 
sampling period 1 was approximately 3 meters per second and gradually 
increased to about 3.5 meters per second at the completion of sampling 
period 3. Hydrostatic conditions were unstable since the water surface 
temperature was 3 2 higher than the air temperature at 4.0 meters. 
During sampling period 4 the wind speed was higher, approximately 
5.5 meters per second at the 4 meter level. Hydrostatic conditions 
were also unstable during this period with the water temperature being 


about 5°C higher thanpthe air temperature at 4. 0 meters. 


fae DISPLAY OF DATA 

CAL-COMP piots were obtained of ail daia sets to facilitate plume 
recognition. Two plots, Seas to the two heights at which the 
Gata were recorded, orem traced on each graph. Each graph was 
fifteen inches in length and contained nine hundred points per plot 
tare represented three minutes of data. Six graphs were needed, 
therefore, to plot data from each eighteen minute period. 

Traces defining the standard plume were then correlated with their 
respective temperature traces in each data set producing a pair of 
correlation curves for each data set. The sum and product of these 
two correlation curves were also obtained. Comparison of all cor- 
relation curves thus obtained showed that the product curve provided 
the best display for plume detection. Typical examples of the tem- 
perature traces investigated and their corresponding product of cor- 


relation curves are reproduced in Figures 12, 13, and 14. 
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Paeure |2. Examples of Traces from Training Set Data. 
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figure 13. Examples of traces obtained from first test set data. 
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Figure 14. Examples of traces obtained from second test set data. 
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C. TABULAMIC Or hike Giles 
Using the scale on the product of correlation curve, the following 


thresholds for various probabilities of plume detection were assigned: 


Low probability - product curve peaks between 
300 units and 400 units, 


Medium probability - product curve peaks between 
400 units and 800 units, 


High probability - product curve peaks above 
800 units. 
The accuracy of each level of detection probability for the training set 
data and the two test sets of data is summarized in Tables III, IV, and 


V respectively. 


oe 





Probability Number of Correct . Incorrect Percentage 


of Comrect Pinas Identifi- Identifi- of 
Dereckron imdicated cation cation Accuracy 
LOW 14 8 6 oy ea 
MEDIUM 9 6 3 66.7 
HIGH ©. 3 0 100. 0 
Totals 26 17 9 65.4 
SUMMARY: 


an 
_— 
wen ao 


There were a total of 23 plumes visually observed 
in this data set of which 17 were detected on the 
correlation plots. The overall percentage of de- 
tection was 74. 0%. 


{Z) The total number of incorrect plume detections 


was. The percentage of incemmect deteetions 
by probability category is as follows: 


Low probability: 66. 7% 
Medium probability: 33.3% 
High probability: 0% 


(3) A total of 6 plumes were visually observed which 
were not detected by the correlation plots. 


Table III. Analysis of Results on Training Set Data. 


Sto 








Probability 
of Correct 
Detection 


LOW 


MEDIUM 


HIGH 


Totals 


SUMN 


(1) 


(3) 


maple IV, 


Niimanper Of iG@ourect inieorne ct Percentage 
Plumes Identifi - Identifi- on 
Indicated cation cation Accuracy 
26 12 14 46. 4 
ad Li 10 Sin 18) 
4 os 0 ICL) Opal 6 
2 f/ ae 24 Oo ae 


There were a total of 45 plumes visually observed 
in this data set of which 33 were detected on the 
correlation plots. The overall percentage of de- 
tection was 73.3%. 


The total number of incorrect plume detections 
was 24. The percentage of incorrect detections 
by probability category is as follows: 

Low probability: 58.3% 

Medium probability: 41.7% 

High probability: 0% 


A total of 12 plumes were visually observed which 
Memcenoe Cerecteauny tne corre raion: plots. 


Analysis of Results on First Test Set Data. 
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Probability 
of Correct 
Detection 


LOW 


MEDIUM 


HIGH 


Totals 


SUMMARY: 


(1) 


(3) 


Nia Deiat NOt neces Incorrect Percentage 


Plumes Identifi- Identifi- Ot 
Indicated cation cation ACCMIracy 

Zn) 1] 10 SS r 

(315; 18 i 15) 

Z Z, 0 100. 0 

48 | eal Oe 64.5 


There were a total of 42 plumes visualiy observed 
in this data set of which 31 were detected on the 
correlation plots. The overall percentage of de- 
tection was 73.1%. 


iheptotalmumber of incorrect plume devections 
was 17. The percentage of incorrect detections 
by probability category is as follows: 

Low probability: 58. 8% 

Medium probability: 41.2% 

High probability: 0% 


A total of 11 plumes were visually observed which 
were not detected by the correlation plots. 


Table V. Analysis of Results on Second Test Set Data. 
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VI, CONCLUSIONS 


The percentage of ome plume detections provides justification 
for the selected correlation approach. Accuracy could be increased 
by further refinement of the standard plume. This could be accom- 
plished by using more than ten known plumes to form the standard. 
Further extensions of this research might be to examine the possibility 
of using other pattern recognition techniques for improved plume 
detection. As an example, an adaptive approach could be employed 
whereby a standard plume would be modified by each newly-identified 
plume to provide the most representative or current pattern. Sucha 
scheme might allow for plume detection under different weather con- 

. ditions. However, the methods and results developed in this study 
appear to be quite useful in their present form. 

Most procedures now used to study atmospheric turbulence depend 
upon assumptions that the processes are nearly random and stationary. 
These assumptions often lead to the use of spectral methods for ex- 
amining the relationship between parameters contributing to heat flux 
(i.e., wand T) or momentum flux (i.e., uandw). The data being 
considered herein, however, show that a significant and persistent 
mecurrence in supposedly random temperature fluctuations has been 
organized plumes. These are intermittent in their occurrence yet all 
have similar features which are identifiable. The intermittency of 
these plumes precludes the use of most techniques employed in atmos- 
pheric turbulence investigations for their study. In order to more 


accurately define the processes responsible for a significant amount 
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of heat flux near the surface, inthis case over a water surface, new 
techniques suchas pattern recognition methods are needed. 

In conclusion, this study has yielded a new ''coordinate" system 
with which to examine processes contributing to the exchange of 


sensible heat between the ocean and the atmosphere. 
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